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ABSTRACT 

The observed masses of the most massive stars do not surpass about 150 Mq. This 
may either be a fundamental upper mass limit which is defined by the physics of 
massive stars and/or their formation, or it may simply reflect the increasing sparsity 
of such very massive stars so that observing even higher-mass stars becomes unlikely 
in the Galaxy and the Magellanic Clouds. It is shown here that if the stellar initial 
mass function (IMF) is a power-law with a Salpeter exponent (a — 2.35) for massive 
stars then the richest very young cluster R136 seen in the Large Magellanic Cloud 
(LMC) should contain stars with masses larger than 750 Mq. If, however, the IMF 
is formulated by consistently incorporating a fundamental upper mass limit then the 
observed upper mass limit is arrived at readily even if the IMF is invariant. An explicit 
turn-down or cutoff of the IMF near 150 Mq is not required; our formulation of the 
problem contains this implicitly. We are therefore led to conclude that a fundamental 
maximum stellar mass near 150 Mq exists, unless the true IMF has a > 2.8. 
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1 INTRODUCTION 

The question on the existence of a finite stellar upper 
mass limit has a long history of debate in the literature 
femiegreenll200(]| ; IMassevI Il99ct references therein). Obser- 
vational evidence for such a limit is scarce because stars 
more massive than 60 — 80 Mq are very rare. While stel- 
lar formation models lead to a mass limit near 100 Mq im- 
posed by feedback on a spherical accr etion envelope jKahnl 
Il974l: IWolfire fc CassmelTilll986l . Il987ft . theoretical work on 
the formation of massive stars through disk-accretion with 
high accretion rates th ereby allowing thermal radiation to 
escape pole-wards (e.g.. lNakanoll989l : IJiiina fc AdamJl996ll 
call the existence of such a limit into question. Some mas- 
sive stars may also form by coagulation of intermediate- 
mass proto-stars in very dense cores of emerging embed- 
ded clusters driven by core-contraction due to very rapid 
accretion of gas with low specific angular momentum, 
thus again avoiding the theoretical feedback-induced mass 
limit dBonnell. Bate fc Zinneckeilll998t Tstahler. Palla fc Hoi 
I2000D . 

In his review, iMassevI (I1998T) points out that inferring 
the masses of very massive stars is difficult due to the fact 



* E-mail: weidner@astrophysik.uni-kiel.de 
f E-mail: pavel@astrophysik.uni-kiel.de 
| Heisenberg Fellow 



that stars heavier than 100 Mq do not have their maximum 
luminosity in the optical bands and are therefore not easily 
discriminated on the basis of photometry from stars with 
somewhat lower masse s. Using combined photo metric and 
spectroscopic methods, Massev fc Hunter! IJ1998) find stars 
with masses ranging up to m = 140 Mq (or even 155 Mq 
depending on the stellar models used) in the rich (about 
10 5 stars) and very young (1-3 Myr) R136 cluster in the 
Large Magellanic Cloud, and that the IMF has a Salpeter 
exponent (a = 2.35) for 3<m/M e < 100. Given this IMF, 
IMassevI il99STl emphasises that the observed most-massive- 
star-mass of around 150 Mq is simply a result of the extreme 
rarity of even more massive stars, rather than reflecting a 
fundamental maximum stellar mass: the observed numbers 
of very massive stars are consistent with the numbers ex- 
pected from sampling from the IMF and the number of stars 
in a cluster. 

In order to re-address this last point, we t ake a n ap- 
proach similar to the route taken bv lElmegreenl <l2000l) . but 
we rely on a different mathematical formulation. The idea 
is to quantify the expected mass of the most massive star, 
m m „, as a function of the stellar mass, M ec i, in an embed- 
ded cluster, and to show that very rich clusters would pre- 
dict an m max which is significantly larger than the observed 
most massive star. Thus we adopt the observed IMF and 
demonstrate that the observed cutoff mass is significantly 
below the expected maximum stellar mass in rich clusters 
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Figure 1. The "logarithmic" IMF (i;h(m) = §(m)m In 10 over 
logarithmic stellar mass above 80 Mq for three different cases. 
The solid line shows an unlimited Salpeter IMF, the dotted line a 
Salpeter IMF truncated at 150 Mq and the dashed line a Salpeter 
IMF limited at 150 Mq = m max * in a way described further in 
§ l2l All three cases are normalised to the same area over 0.01 ^ 
tu/Mq < 00. 



if there were no fundamental upper mass limit. The impli- 
cation would thus be that there must exist a fundamental 
upper mass limit, m max «, such that m max < m max * for all 
M ec i. With the use of simple equations concerning the IMF, 
and the real ization that most if not all stars are born in stel- 
lar clusters jLada fc Ladall2003D with an universal IMF, we 
show that the solutions of these equations predict a very dif- 
ferent high mass spectrum for a finite or infinite fundamental 
upper stellar mass, m max *, in dependence of the associated 
cluster mass. The principles are shown in Fig. Q 

The next Section |5] introduces the equations and the 
analytical and numerical methods used to solve them, while 
the results are shown in § [21 The implications are discussed 

in §m 



2 METHOD 

For our calculations we use a 4-component power-law IMF, 



-a 



£(m) 



"Mow ^ m < m H , 
ma sj m < mo , 



( m H ) 



with exponents 

a = +0.30 
ai = +1.30 
a 2 = +2.30 
a 3 = +2.35 



0.01 < m/M Q < 0.08, 
0.08 < m/M Q < 0.50, 
0.50 < m/M Q < 1.00, 
1.00 < m/MQ, 



where dN = £(m)dm is the number of stars in the mass 
interval m to ra + dm. The ex ponents Qj represen t the 
Galactic-field (or standard) IMF ilKroupall200ll. I2002T) . The 
advantage of such a multi-part power-law description are 
the easy integrability and, more importantly, that differ- 
ent parts of the IMF can be changed readily without affect- 
ing other parts. For example, the stellar luminosity func- 
tion for late-type stars p oses significant constraints on 
the IMF below m<lM n jKroupa. Tout fc Gilmord 11993 ; 
iReid. Gizis fc Hawlevl 120021 lKroupall2002ll which therefore 
must remain unaffected when changing the IMF for mas- 
sive stars. The observed IMF is today understood to be an 
invariant Salpeter power-law above a few Mq, being inde- 
pendent of the cluster de nsity and metalicity for metalicities 
Z> 0.002 iMassevlll99sh . 

The basic assumption underlying our approach is the 
notion that stars in every cluster follow this same universal 
IMF. 



2.1 Number of stars 
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The Number of stars above a mass m is 



N 



£(ra) dm, 



(3) 



when the normalisation constant k (eq. is given by the 
stellar mass of the cluster, 



M ficl = 



m • £(m) dm. 



(4) 



Here we use the cluster mass in stars prior to gas-blow-out 
and thus prior to any losses to the stellar population due to 
cluster expansion jKroupa fc Boilvll2002l) . 

In Fig. |21 it is shown that a significant number of stars 
with masses m > 150 Mq should be present in R136 (10 
stars for Mri36 = 5 x 10 4 Mq and 40 stars for Mri36 = 
2.5 x 10 Mq) if no fundamental upper mass limit exists 
(m max * = 00) and if the IMF is a Salpeter power-law above 
about 1 Mq , whereas none are observed. This sets the prob- 
lem for which we seek a solution by considering a finite 
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2.2 The limited case 

First we examine the case were a finite upper mass limit for 
stars exists. Here two upper mass limits have to be differ- 
entiated: the fundamental maximum-mass a star can have 
under any circumstances, m max * , and the 'local' upper mass 
limit m max ^ m max * for stars in a cluster with a stellar mass 
Mod. The mass of the heaviest star in a cluster, m max , fol- 
lows from stating that there is exactly one such star in the 
cluster, 



£(m) dm. 



(5) 



Note that lElmeereenl (I200CT) uses m max * = 00 in his formu- 
lation of the problem. After inserting eq.|5]the integral can 
be solved obtaining (a; 7^ 1): 
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Figure 2. Number of stars (logarithmic) above mass m for R136 
with different mass estimates (dott ed line: A/risr = 2,5 X 10 5 Mq, 
dashed line: A/ H ,i 36 = 5 X 10 4 M^. ISelman et alj|l999l) . The ver- 
tical solid line marks m = 150 Mq. 
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Figure 3. Dependence of the stellar upper mass limit, m ma x, on 
the cluster mass for a limited (?Ti max * = 150 M Q ; dashed line) 
and an unlimited (m max * = oo; dotted line) fundamental upper 
stellar mass and ag = 2.35. 



as long as m max > mi. For mo ^ m max < mi we would 
have 



1= ^u^r^'iC 
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and so on. For the numerical results obtained in this work 
m ma x* = 150 Mq is assumed. 

In order to solve this equation with two unknowns, k 
and m ma x, we need another equation. It is provided by the 
mass in embedded-cluster stars (eq.^J. With the use of £(m) 
(eq.0, eq. El leads to (on / 2): 
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for m max > mi and with mi ow set to 0. 
lows. For mo ^ m max < mi eq. |8]would 
earlier term, and so on. 

Finally inserting eq.||j] after a short 
Ogives M cc i in dependence of '72 max • This 
for m max in dependence of M ec i. This 
the roots of this result after subtracting 
the solution for a power-law with Q3 = 
150 Mq as a dashed line. 



01 Mq in what fol- 
be truncated at an 

transformation into 
must now be solved 
is done by finding 
Med- Fig. El shows 
2.35 and m max * = 



2.3 The unlimited case 

In the case of m max * = 00 eqs 2] and [8] remain as they are 
while and HJ change to 



I = / C( m ) dm 
and (as long as m max > mi and ay > 1) 



1-03 
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respectively. As only the normalisation factor k deduced 
from llOl changes. eq. [HJstays the same, and inserting 1 101 into 
gives M eo i in dependence of m max for the unlimited case. 

Fig. |3 shows that the solution for unlimited stellar 
masses (dotted line) has a much faster rise than the limited 
case. If there were no fundamental upper mass limit for stars 
then a Salpeter IMF would predict stars with much larger 
masses (m max > 200 Mq) for clusters with M oc i > 10 4 ' 5 Mq 
than ar e observed to be p resent. This is also found to be the 
case bv lElmeereenl J2000). 



3 RESULTS 

The results of solving m max (M oc i) for a grid of cluster 
masses ranging from M ec i = 5 Mq (Taurus- Auriga-like stel- 
lar groups) to 10 7 Mq (very massive stellar super clusters) 
are plotted in Figs 2] and 

Fig. 21 shows the variation of the maximum possi- 
ble mass for a star, m max , in dependence of the clus- 
ter mass, M ec i. In the unlimited case (long-dashed line) 
a linear relation (in double logarithmic units) is seen. 
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Figure 4. Double logarithmic plot of the maximal stellar mass 
versus cluster mass. Shown are three cases: finite total upper 
mass limit of m max , = 150 Mq (dotted line), m max * = 1000 Mq 
(short-dashed) and no limit, m max » = oo (long-dashed). The ver- 
tical lines mark the empirical mass interval for R136 in the LMC. 
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Figure 5. Maximal stellar mass versus cluster mass (logarith- 
mic) . Results are shown for different exponents (0:3 ) above 1 Mq 
and for the limited (m max * = 150 Mq) and unlimited case. The 
vertical lines mark the empirical mass interval for R136 in the 
LMC. 



Two vertical lines indicate the observation al mass inter- 
val fo r R136 in the Large Magellanic Cloud llSelman et alJ 
1999). Without a fundamental upper mass limit R136, for 
which iMassev fc Hunter! l|l99£l ) measure a Salpeter power- 
law IMF for m > few Mq , should have stars with m > 
750 Mq, whereas no stars with m > 150 Mq are seen. Sim- 
ilar values are fo und from statistical sampling of the IMF 
felmegreenll2000l) . For m max * = 150 Mq (dotted line), on 
the other hand, the cluster has an upper limit of 140 — 
150 Mq, in agreement with the observational limit. 

The influence of the high-mass exponent 0:3 on the 
m max (M ec i) relation is shown in Fig. [K] Plotted are graphs 
for limited (150 Mq) and unlimited cases, each for 03 = 2.35 
(Salpeter), 2.70 and 3.00. Exponents az > 2.8 lead to a 
w max (M ec i) relation which allows upper masses in R136 of 
around 150 Mq even for the unlimited case (m max * = 00). 
Fig. HJ shows that in the case of R136 and for 03 > 2.8 no 
distinction can be made between m max , = 150 Mq and 00 
given the uncertainty in M ec i. 

Because massive stars are very rare the IMF exponent 
is often based on limited statistics and usually only for stars 
with m 2s 40 Mq. We therefore also consider now the possi- 
bility that the IMF slope is Salpeter to a certain limit (e.g. 
40 Mq) but then turns down sharply. For this purpose we 
set mi = murder in eq.^with Q2 = 2.35 (0.5 Mq - rrtbordor) 
and find that a m>mbord(!r = Q3 such that eq. |5]is fulfilled for 
w max = 150 Mq. The result is plotted in Fig. Q 

From Fig. [T] it is evident that in order to reproduce the 
observed limit of about 150 Mq for R136 from a formally 
unlimited mass-scale and a down-turn mass (mb or dcr) of, 
say, 40 Mq the exponent has to change to a TO >m bordel =3.6 
(for Mri36 = 5 x 10 4 M©) or 4.5 (M R i 36 = 2.5 x°10 5 Mq). 



Such a down-turn near 40 Mq is not seen in those popula- 
tions that do contai n more-massive stars (e.g . R136 contains 
about 40 03 stars. Ma ssev fc Hunterlll998T) . and we there- 
fore consider m max * w 150 Mq as being the more realistic 
possibility. Note though that the existence of m max „ leads 
to a sharp decline of the IMF near 120 Mq which leads to 
a similar effect as an increase of a m > mbordel near this mass 
(Fig-Q- However, our formulation needs one additional pa- 
rameter (m max ») to implicitly account for this turn-down of 
the IMF, while modelling an explicit turn-down would need 
two additional parameters (mborder and Qm>m bordor ). 

For massive stars the multiplicity proportion is typically 
very hig h with most Q st a rs having mor e than one compan- 
ion (e.g. IZinneckeil 120031 : iKroupal 120031) possibly implying 
the true underlying bin ary-corrected IMF to have «3 ^ 2.7 
JSagar fc Richtlerll99ll Weidner & Kroupa, in preparation) . 
If this is the case then m max * cannot be constrained given 
the available stellar samples because the Local Group does 
not contain sufficiently massive, young clusters. 



4 DISCUSSION AND CONCLUSIONS 

With a rather simple formalism based on the current knowl- 
edge of the IMF we have shown that the mere existence 
of a fundamental upper mass limit implies the highest 
mass a star can have in a massive cluster to be differ- 
ent to the case without such a limit. For low-mass clusters 
(Med < 10 3 Mq) the differences of the solutions are neg- 
ligible (Fig. |1J, but in the regime of the so-called 'stellar 
super-clusters' (M ec i > 10 4 Mq) they become very large. 
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Figure 6. The mass limits (m max ) in dependence of the IMF 
exponent «3 (above 1 Mq) in the limited case (m max » = 150 Mq) 
and the unlimited case (m max * = oo) for the two mass limits of 
R136 shown in Figs. Ill and l5l 
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Figure 7. The power-law exponent a needed to produce a high- 
mass limit of 150 M© for R 136 (solid line: M R136 = 2.5 X 10 5 Mq 
and dotted line: M R i 36 = 5 X 10 4 Mq) when the IMF is Salpeter 
up to a certain mass limit mborder- 



Without such a limit, clusters like R136 in the LMC would 
hav e stars with m > 7 50 Mq . 

lElmegreenl (|200(t l presents a random sampling model 
for star formation from the IMF which is similar to our 
model. However, Elmegreen assumes a Salpeter power-law 
IMF above 0.5 Mq and no specific stellar mass limit. In order 
to reduce the number of high-mass stars above ~ 130 Mq he 
assumes an exponential decline for the probability to form 
a star afte r a t urbulent crossing time. The results of the 
lElmegreenl i2000T) model are summarised by him as follows: 
"There is a problem getting both the Salpeter function out 
to ~ 130 Mq in dense clusters and at the same time not 
getting any ~ 300 Mq stars at all in a whole galaxy." 

He discusses the following six explanations for this prob- 
lem: 

i. Stars more massive than ~ 150 Mq exist but have not 
been found yet. 

ii. A self-limitation in the star formation process prohibits 
stars above a certain limit. 

iii. Super-massive stars exist but evolve so quickly that 
they do not leave their primordial clouds - making them 
observable only as ultra-luminous infrared sources. 

iv. An assumed limit of the cloud size for coherent star 
formation 

v. The star forming clouds are destroyed after a star of a 
certain (maximum) mass forms. 

vi. The IMF is not universal but different for various star 
forming regions. 

Case i can be excluded here because of the number of 
super-massive stars expected, for example in R136. Concern- 
ing case iii no such sources have been found to our knowl- 
edge. The cases ii, iv and v lead to a physical upper limit 
consistent with this work. From the point of view of this 
work it is not possible to differentiate between them. Fi- 
nally as several observations of various cl usters show a uni- 
versa l Salpeter IMF up to ~ 120 Mq (e.g.lMassev fc Hunterl 
Il99ct ISelman et al.lll999t ISmith fc Garlagherll200ll) case vi 
appears unlikely. Elmegreen thus sees the finite upper mass 
limit as a cut-off to the unlimited solution. 

In contrast, we introduce the fundamental upper mass 
limit consistently into the formulation of the problem, and 
together with the use of a realistic IMF we are able to show 
strong deviations of the solutions beyond a simple cut-off. 
The formulation presented here has the advantage of ex- 
plaining the observations under the rather simple notion 
that all stars form with the same universal IMF. 
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